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Abstract. The Auger Collaboration reports 1, 2J that the arrival directions of ^ 60 
EeV ultra-high energy cosmic rays (UHECRs) cluster along the supergalactic plane 
and correlate with active galactic nuclei (AGN) within « 100 Mpc. The association of 
several events with the nearby radio galaxy Centaurus A supports the paradigm that 
UHECRs are powered by supermassive black-hole engines and accelerated to ultra- 
high energies in the shocks formed by variable plasma winds in the inner jets of radio 
galaxies. The GZK horizon length of 75 EeV UHECR protons is « 100 Mpc, so that 
the Auger results are consistent with an assumed proton composition of the UHECRs. 
In this scenario, the sources of UHECRs are FR H radio galaxies and FR I galaxies like 
Cen A with scattered radiation fields that enhance UHECR neutral-beam production. 
Radio galaxies with jets pointed away from us can still be observed as UHECR sources 
due to deflection of UHECRs by magnetic fields in the radio lobes of these galaxies. 
A broadband ~ 1 MeV - 10 EeV radiation component in the spectra of blazar AGN 
is formed by UHECR-induced cascade radiation in the extragalactic background light 
(EEL). This emission is too faint to be seen from Cen A, but could be detected from 
more luminous blazars. 



PACS numbers: 98.54.-h, 98.62.Js, 98.70.Sa 
Keywords: Ultra-high energy cosmic rays; radio galaxies; active galactic nuclei, black 
holes 

Submitted to: New J. Phys. 



Ultra high energy cosmic rays from black hole jets of radio galaxies 



2 



1. Introduction 

The Auger Observatory in tlie Mendoza Province of Argentina at ~ 36° S latitude 
determines the arrival directions and energies of UHECRs using four telescope arrays to 
measure Ni air fluorescence and 1600 surface detectors spaced 1.5 km apart to measure 
muons formed in cosmic-ray induced showers. Event reconstruction using the hybrid 
technique gives arrival directions better than 1°, and energy uncertainties at 10^° eV (100 
EeV) of ~ 11% for a 50% Fe and 50% p composition [3]. Analysis of the composition 
of the high energy showers in the early Auger analysis showed it becoming heavier, 
somewhere between p and Fe, at lO^^'"^ eV |ll[5]. By contrast, HiRes data are consistent 
with dominant proton composition at these energies [6] , but uncertainties in the shower 
properties [7] and particle physics extrapolated to this extreme energy scale [8] preclude 
definite statements about composition. 

In the Auger analysis [H [2], a probability statistic P corrected for exposure is 
constructed from the nearest-neighbor angular separation ip between the arrival direction 
of an UHECR with energy E and the directions to AGN in the Veron-Cetty and Veron 
(VCV) catalog [9], containing 694 active galaxies with z < 0.024 or distance d < 100 
Mpc. P was minimized for ip = 3.1°, threshold clustering energy Ed = 56 EeV, and 
clustering redshift Zd = 0.018 {d^ = 75 Mpc), containing 27 events (the two highest 
energy events were 90 and 148 EeV). Twelve events correlate within 3.1° of the selected 
d < 75 Mpc AGN, and another three within the vicinity of one of these nearby AGN, 
ruling out an isotropic UHECR flux or a Galactic source population. Note that the 
VCV AGNs are not necessarily the sources of the UHECRs, but may only trace the 
same matter distribution as the actual UHECR sources. 

This discovery opens the field of charged-particle astronomy. Only at the highest 
energies can arrival directions of charged particles be associated with their sources. 
This is because deflections by (i) the Galactic or (ii) intergalactic magnetic (IGM) field 
isotropizes the directions of lower energy cosmic rays. 

Here we consider the Auger clustering results within the paradigm that extragalactic 
black-hole jet sources accelerate UHECRs. Estimates of deflection angle and delays are 
given in Section 2, specialized to an assumed proton composition of the UHECRs. The 
question of the GZK cutoff and the UHECR horizon is revisited in Section 3. In Section 
4, we apply standard synchrotron theory to the lobes of Cen A in order to estimate the 
equipartition magnetic field and absolute jet power, using a new technique that uses 
the jet /counter-jet ratio to determine the speed of the outflow. From this, limits on 
UHECR acceleration in colliding shells of blazars are used to derive maximum particle 
energies. Also, an estimate is made of the flux of Compton-scattered cosmic microwave 
(blackbody) background radiation (CMBR), and compared with the flux of secondary 
nuclear production in Cen A's radio lobes. In Section 5, we present calculations of the 
broadband 7-ray uFu flux from Cen A due to secondary cascading of protons on the 
EBL, for 1 nC (= 10"^ G) IGM fields and show that Cen A is not detectable with 
current instrumentation. Concluding remarks are given in Section 6. 
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2. Magnetic Field Deflections of UHECRs 

For case (i), the Galactic magnetic field can be approximated by a magnetic disk with 
characteristic height hmd, giving a deflection angle 9dfi ~ hmdcscb/ri^. Here b is the 
Galactic latitude of the UHECR source, and the Larmor radius of a particle with energy 
E and charge Ze is 

(E/60EeV) , 

Thus 

ZhrrJ^pc)B(jlG) 

^'^f^'^^^-^ sin6(E/60EeV) ' 

limited by the finite extent of the magnetic disk. Mean magnetic fields B in the ^ 0.2 
kpc thick gaseous disk of the Galaxy are ~ 3 - 5 fiG, but could fall to <^ lyuG in the 
kpc-scale halo [10]. Deflection angles ^ 3° from Gen A (6 = 19.4°, galactic longitude 
£ = 309.5°, declination —43°, distance d = 3.5 Mpc) restrict UHEGRs to protons or 
light-Z nuclei and a small ( ^ 0.1 fiG) Galactic halo magnetic field. 

For case (ii), the defiection angle of an UHEGR ion when propagating through 
the IGM field from a source at distance d is [HI [12] OdjcM — djlri^^fW^^ ~ 
0.04° Z (5_i2) rf(100 Mpc)/[E(60 EeV)yi\Q , where A^i„^ ~ max(d/A, 1) is the 
number of reversals of the magnetic field (also expressed through the magnetic-field 
correlation length A) and (-B-12) is the mean magnetic field of the IGM in pico-Gauss 
(1 pG = 10~^^ G). If UHEGRs within 3° of Gen A are accelerated by the radio jets of 
Gentaurus A, then to avoid much larger defiections than those made by the Galactic 
magnetic field requires that (5_i2) < 2000 ^fN~ E(60 EeV)/[Zrf(3.5 Mpc)]. By this 
reasoning, the mean IGM field in the directions towards AGNs 75 Mpc distant in the 
supergalactic plane (SGP) is restricted to be ^ lOOv^A^ i?(60 EeV)/Z pG [13j. Time 
delays between electromagnetic outbursts and UHEGR arrival windows from Gen A due 
to propagation through the IGM are about 



^ n ^ ^^(3-5 Mpc)Z2(5_i2) 
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At ^ ^ 0.5 —7- ^ , day , (3) 

where the final expression holds because Nmv ~ 1- Variable 7-ray fiaring activity from 
Gen A could be refiected in variable UHEGR activity on sub-day timescales for IGM 
fields (-B-12) ~ 10 and A 0.1 Mpc or a large-scale (~ Mpc) ordered field with mean 
strength (-8^12) ~ 1- 

The most intense magnetic fields between us and Gen A consistent with the 
defiection data for an assumed proton composition of the UHEGRs are {B) ^ 2 y/Ninv 
nG. For the Auger data reaching to galaxies at distance ci ^ 75 Mpc, {B) ^ 100 y/N^ 
pG. For Ninv ~ 100 (A = 1 Mpc), (B) ^ 1 nG. For the calculations in Section 5, we use 
(B) = 1 nG which is a more realistic value because UHEGRs are correlated with more 
distant sources as well. 



Ultra high energy cosmic rays from black hole jets of radio galaxies 



o 

Q. 



(/} 
(O 
O 
_l 

0) 

c 
LU 




Figure 1. Mean- free paths 
for energy loss of UHECR 
protons in different model 
EBLs are shown by the 
solid curves, with pho- 
topair (dotted) and photo- 
pion (dashed) components 
shown separately. "CMB 
only" refers to total energy 
losses with CMB photons 
only. Inset: Measurements 
of the EBL at optical and 
infrared frequencies, includ- 
ing phenomenological fits to 
low-redshift EBL in terms 
of a superposition of mod- 
ified blackbodies. A Hub- 
ble constant of 72 km s^^ 
Mpc~^ is used throughout. 
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3. UHECR Proton Horizon 

The clustering energy Ed = 60 EeV separates UHECRs formed mainly by sources along 
the SGP at c? dd from lower-energy UHECRs formed nearby and on ^ 75 Mpc scales. 
A high-significance steepening in the UHECR spectrum ai E = 10^^'^ eV = 4 x 10^^ 
eV and at E = 10^^-^ eV = 6 x 10^^ eV was reported, respectively, by the Auger [4] 
and HiRes [14] collaborations in 2007. These results confirm the prediction of Greisen, 
Zatsepin and Kuzmin fTEl [TB] that interactions of UHECRs with CMBR photons cause 
a break in the UHECR spectral intensity near 10^° eV. 

The clustering and GZK energies coincide because ^ 60 EeV UHECR protons 
originating from sources at the ^ 100 Mpc scale have lost a significant fraction of 
their energy due to photopion losses with the CMBR, so that higher-energy particles 
from the more distant universe cannot reach us. The energy-loss mean free path 
r^^ = ct^T^ = C'y\d'y / dt\'^l of an UHECR proton with energy E = 10^°i?2o eV to 
photopion losses with the CMBR at low redshifts is given, in good agreement with 
numerical calculations [HI [IE] , by the expression 

, 13.7exp[4/E2ol 
^*-<^-' ^ 11 + 4/4] """^ 
. The term r^^ gives the mean distance over which a particle with energy E loses 
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^ 1 — = 63% of its energy. The MFPs for energy loss by photopair and photopion 
losses in different model EBLs, including an EBL consisting of the CMB alone, is shown 
in Fig. [H 




Energy (eV) 

Figure 2. Heavy dotted curves give tlie iiorizon distance for UHECR protons as 
a function of total proton energy, using the local CMBR and the low and high EBL 
target radiation field shown in the inset to Fig.[TJ Analytic approximation to photopion 
energy loss mean- free path, eq. is given by dot-dashed curve, and the UHECR 
proton horizon, eq. ([5]), by dashed curve, for CMBR only. Auger clustering results 
are indicated by the shaded arrows. Distances to key radio galaxies are shown, and 
light dashed lines separate quasi-linear proton trajectories from trajectories with strong 
deflections in overall ordered magnetic fields with strength as labeled. Short-dashed 
lines give length for a proton to be deflected by 0.16'_i rad in magnetic fields ranging 
from fO-s - 10-" G. 

The GZK horizon length giving the mean distance from which protons detected 
with energy 10^°£'2o eV originate depends in general on injection spectra and source 
evolution [19] (see also [20l [21]), but a model-independent definition that reduces to 
r^n{E2o) for an energy-independent energy- loss rate considers the average distance from 
which a proton with measured energy E had energy cE. The horizon distance, defined 
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this way, is given by 

rhrz{E2o) = / — r^^{x) = 

JE20 X 

13 7 dx "^P^^/^^ Mdc = 1-1^20 ^^P(4/^^o) Mnc (5) 
Je2o x{1 + 4/x) ^ 1 + 1.6E|o/13.7 ^ ^ 

where the last two expressions give the proton horizon on CMBR photons alone. 

Using the phenomenological fits to low and high EBLs at optical and IR frequencies, 
represented in the inset to Fig. [T] as a superposition of blackbodies [T3], gives the 
corrected horizon distance shown in Fig. [2l The horizon distance for ^ 57(75) EeV 
UHECR protons is ^ 200(100) Mpc. A GZK horizon smaller than 40 Mpc applies to 
protons with E ^ 100 EeV. This explains the clustering observations observed by the 
Auger collaboration if UHECRs are predominantly protons, but would be inconsistent 
with the Auger results if the UHECRs are composed of high-Z material, as the mean free 
path for photodisintegration can be much larger than 100 Mpc for ^ 10^° eV ions like 
Fe [13]. If the Auger energy scale is underestimated by ^ 20% due to systematic effects 
(which would also reconcile the discrepancy between the different HiRes and Auger GZK 
energies as determined from the spectral break), then GZK losses on UHECRs with a 
dominant proton composition would explain the clustering towards the SGP even more 
decisively. The higher energy scale for the Auger experiment might also explain the 
lack of clustering observed by HiRes at > 57 EeV [22], which would at this energy cut 
in the HiRes experiment include large numbers of lower-energy, more distant and less 
clustered cosmic rays. 



4. UHECRs from AGN Jets 



Discovery of UHECR arrival directions clustered towards the SGP was anticipated by 
analysis of UHECR data from Haverah Park, AGASA, Volcano Ranch, and Yakutsk 
observatories [23]; see also [ITJ. Compared to an isotropic source flux, the average and 
rms angular distances toward the SGP were enhanced at the 2.5 - 2.8 a level for events 
with E > 40 EeV. Stanev et al. [23] argued that their analysis favors radio galaxies as 
the sources of UHECRs. A radio-galaxy origin of UHECRs is consistent with harder 
radio sources [21] and clustering of sources in the two Jansky, 2.7 GHz Wall and Peacock 
catalog [25] towards the SGP. 

The SGP runs through the Virgo Cluster at ^ 20 Mpc, and contains an assortment 
of radio galaxies such as M87, Cen A and NGC 315, and the starburst galaxies M82 
and NGC 253. Using infrared galaxy surveys to better define the SGP improves the 
significance of correlations of UHECRs with the SGP [26]. When weighted by hard 
X-ray flux, the UHECR arrival directions are strongly correlated with Swift Burst Alert 
Telescope galaxies within 100 Mpc, which trace the SGP [27] . Searches [28l [29] for 
specific AGN in the VCV and NASA/IPAC NED catalogs finds Seyfert 2, low ionization, 
and other radio-quiet galaxies closest to the UHECR arrival directions, in addition to 
associations with Cen A, Cen B, an FR II radio galaxy and a BL Lac object within 140 
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Mpc. IGR J21247+5058, an FR II broad-lined radio galaxy at z = 0.02 or ^ 80 Mpc, 
recently discovered with INTEGRAL [30], is 2.1 degrees away from a HiRes Stereo event 
with ^ > 56 EeV [31]. At least 8 of the 27 UHECRs with E > 56 EeV are within 3.5° 
of nearby radio galaxies [32] . 

We suppose that the evidence is compatible with an AGN origin of UHECRs 
in the black-hole jets of radio galaxies, with Centaurus A being the most prominent 
example. Gen A itself, though classified as an FR I radio galaxy, has a bolometric radio 
luminosity exceeding 4 x 10^^ ergs s~^ [33], near the dividing line between FR I and FR 
II galaxies in terms of radio power. The giant elliptical galaxy in Gen A is intercepted 
by a small spiral galaxy making the prominent dust lanes [M] which, including inner 
dust torus emission, could conceal the optical/UV line strengths and distribution of 
broad line region gas near the Gen A nucleus [35l [36] . A strong scattered radiation field 
is important for photomeson losses and neutral beam production in AGN jets [371 EH], 
and seems required in many BL Lac objects because of spectral fitting difficulties with 
synchrotron/SSG models [39lll0lliT]. 

The nucleus of the Gen A jet is visible at radio and X-ray energies, revealing sub- 
luminal {v ~ c/2) relativistic outfiows and jet/counterjet fiuxes of X-ray knots consistent 
with mildly relativistic speeds on projected ~ 10 pc scales [121 113]. Gen A was detected 
at hard X-ray and 7-ray energies with OSSE, GOMPTEL and EGRET on the Compton 
Observatory [^ H5] . emitting ^ 5 x 10^^ ergs s~^ in keV - MeV radiation with day- 
scale variability observed at ~ 100 keV (see also [16]). Most of this emission is probably 
quasi-isotropic radiation from the hot accretion plasma. Variability of ~ GeV radiation 
could be detected from Gen A with the Fermi Gamma ray Space Telescope if the jetted 
7-ray emission is not too dim. Such a detection would also discriminate between an 
inner jet and extended origin of the 7-ray emission, for example, from Gen A's lobes 

4-1. Jet Power from Synchrotron Theory 

The corrected Auger point source exposure is x^{^s) l^m — (9000 x 0.64/7r) km^ yr, 
where x is the exposure, u!{6s) = 0.64 is an exposure correction factor for the declination 
of Gen A, and Qqq = vt is the Auger acceptance solid angle [IE]. For a power law proton 
injection spectrum with number index a, the apparent isotropic UHEGR luminosity 
from Gen A is 

LcenA= 1.6x10^x60°-^^^^^^ (^) [E(EeV)]^-<^ crgs s-^ (6) 

which for = 2 events gives LcenA(a = 2.2) ~ 5 x 10^^[-E'(EeV)]-°-^ ergs s"^ and 
LcenA{a = 2.7) ~ 2 X 10^°[E(EeV)]-°-^ ergs s"! J The Gen A radio luminosity thus 
greatly exceeds the UHEGR luminosity unless cosmic rays are strongly beamed along 

I This power should be increased by factors 2 - 3 if 5 or 6 events are associated with Cen A and its 
lobes [ig. 
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the Cen A jet. The integrated cosmic ray power for a = 2.0 from 1 GeV to 10^*^ eV is 
^ 2 - 3 X 10"^° ergs s~^ for = 5 events. 

It is interesting to note that the production of ~ 10^° ergs s^^ in UHECRs by Cen 
A, the dominant radio galaxy within ~ 10 Mpc, represents an UHECR emissivity of 
~ 8 X 10^3 ergs yr'^ Mpc"^ (cf. [i9]). 

4-1.1. Equipartition Magnetic Field The UHECR luminosity can be compared with 
the time-averaged jet power inferred from the magnetic field energy in the radio lobes, 
assumed to be inflated by the pressure of the black-hole jet. Knowing the power of a 
black- hole jet allows one to derive the maximum particle energy of UHECRs accelerated 
through Fermi processes. 

Here we state some elementary results from synchrotron theory needed to derive 
the jet power [50]. We treat a standard blob model and assume that the measured 
radio emission is nonthermal synchrotron radiation from randomly oriented electrons in 
a randomly directed magnetic field. In this approximation, the equipartition magnetic 
field is defined by equating the magnetic-field energy density Ub = -B^/Stt with the 
total particle energy density Upar consisting of hadrons and electrons, assumed to be 
dominated by a large-scale, randomly oriented magnetic field of strength B. 

The equipartition magnetic field B^q = B{keq = 1) with Upar = keqUs, where 

2/7 



Bik. 



B. 



cr 



eq) 



6d 



37rdimecV.rin(e2/ei)(l + Cpe; 



(7) 



V^caTUlk,,^il + z)e. 

Here B^r = mlc^/eh = 4.414 x 10^^ jj^^ ^ B^/Sn = 7.752 x 10^^ ergs cm'^, and 
Sd is the Doppler factor. In this expression, a factor (p^ more energy is assumed to be 
carried by protons and ions than leptons. This result applies to the Fi, oc z/~^/^ portion 
of the synchrotron spectrum made by electrons with N'l'y') oc 7'"^ that carry most of 
the energy. 

The uF^ synchrotron flux //^" at photon frequency u = meC^e/h is found in eq. 
([7]). Note that f^^"'/y/e2 is constant throughout the range ei < e < £2 when f^^^ oc e^^^, 
so that the value of f^^^/y/e is independent of e over this range. Also, Beq depends 
only very weakly on the normalization of the energy contained in electrons, varying 
oc [ln(e2/ei)]2A. 

4.1.2. Synchrotron Power The total jet power of a one-sided jet, referred to the 
stationary black-hole reference frame, is then 



P* ^ Tir'^pcT^ 



_87r ^ VI 2caTM V £2(1 + 2) ""^ei' 
defining 63 = B / B^r From this, one can show [501 IH] that 

p;{B) = \p;{B„,.nL){u' + ^), (9) 



(8) 



Ultra high energy cosmic rays from black hole jets of radio galaxies 



9 



and u = B/BminL ■ The minimum power is 

P;iB^^nL) = lnc/3T'r','U^r . (10) 

The magnetic field giving the minimum jet power required to account for the synchrotron 
fiux is 

BminL = i^fBe, = OmiBeg . (11) 

4.1.3. Application to Centaurus A The distance to Cen A is dcen a = S.Siis.s Mpc, 
with 6/3.5 — 1- The preceding relations give 

r n 2/7 

PnB„..C, ^ ^..^(f 274m.cVr-.fe/^.)(l + C.) ,,,, 



This applies to either spherical or cubical volumes by writing the emitting volume in 
the fluid frame as 

K = = 6.7 X lo69f7(rf3.5V'o)' cm3 . (13) 

The geometry factor is given hy g = 1 for cubical and g = t^/Q for spherical volumes. 
The angle ipQ, in degrees, is the angular extent of the emitting volume. 



Writing f^"^^ = fe2\J^/^2^ with /^^ = 10~^^/_i2 ergs cm~^ s~^, and normalizing to 
the 21 cm (1.43 GHz) fiux, so €2 = 1.15 x 10^^^e2i, then the magnetic field associated 
with the minimum jet power is 



The minimum jet power for a one-sided jet is, therefore, 

r.2 



P*iB^^nL)=^y<W'''|3[ — ] X 



5d 



(\ 4/7 
tofeAO(i + C) ^^g^/^, (15, 
\/f2l(l + j) / 

Taking the bolometric factor ln(e2/ei)(l + Cpe) = 10 in eqs. ( fT4l) and ( fT5l) gives the 
results shown in Table 1 for the five different regions defined in Ref. [17]. Uncertainties 
by a factor 10 in the bolometric factor translate to a factor 10^/^ = 1.93 in BminL and 
a factor 10^/^ ^ 3.73 in P„anL- 

To get the speed of the outflow, we use the jet/counterjet ratio p, which is a number 
taken directly from observations. If 9 is the angle of the jet nearest to the line of sight, 
and we assume that all jets are two-sided with equal power ejected in the opposite 
direction, then it is easy to show [5ll [52] that the speed 

COS 6 p^i ' + 1 
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Table 1. Jet power of Con A from synchrotron theory** 



Regions'^ 


F(Jy) 


/-12 




g 




^minL 


p;(B,„.„i,) 

(r/5D)2(A/10)4/7 

10^*3 ergs s-i 


1 


91.2± 2.1 


1.3 


2.0 


1 


0.1 


0.97 


2.1 


2 


74.6 ±3.0 


1.5 


1.0 


1 


0.2 


0.99 


3.4 


3 


545 ± 1.3 


7.78 


2.0 


7r/6 




1.96/r 




4 


204.7 ±4.3 


2.92 


1.7 


1 


0.2 


1.26 


6.2 


5 


111.2 ±6 


1.59 


1.6 


1 


0.1 


1.08 


2.2 


Total 


1026.7 ±8.3 


14.6 













^Radio observations at 21 cm from Ref. [47] . 
^Regions defined in Ref. J7J- 



Thus we interpret the asymmetrical lobe flux in Cen A as a consequence of aberration of 
the mildly relativistic two-sided outflow, rather than as a difference of spectral properties 
due to distinct environmental effects. The jet power obtained from this interpretation 
can be compared with other methods to determine jet power [53J. 

Table [1] shows that the jet/counterjet ratio p ^ 1.5 - 2 for regions 1 and 5, and 
p 2 - 4 for regions 2 and 4, each pair being at equal angular separation from Cen 
A's nucleus. The orientation 6 of the jet of Cen A with respect to our line of sight is 
35° ^ 9 ^ 72° [51], implying that /3 ~ 0.1 for regions 1 and 5, and /? ~ 0.2 for regions 
2 and 3, which are the values used in the calculations of the jet power shown in Table 
1. These speeds refer to the flow of plasma on the scale of hundreds of kpc from the 
nucleus of Cen A, so that the plasma from inner jet has to be ejected with much higher, 
probably relativistic speeds. The factor (T/6j))^, when included, reduces the power of 
the brighter lobe and increases the power of the dimmer lobe to reproduce the underlying 
assumption in the method that the powers of the two jets are equal. The minimum total 
jet power to produce the weakly boosted Cen A radio emission is therefore «i 2 x 4 x 10"^^ 
ergs s~^, or total jet power P*°* ^ lO''^ ergs s~^, as follows from Table 1. 

With total mean absolute Cen A jet powers ~ 10^^ ergs s~^, apparent isotropic 
jet powers can reach 10^^ - 10"^^ ergs s~^ during flaring intervals. Indeed, apparent 
isotropic flaring luminosities in 7 rays alone exceed 10^^ ergs s~^ in Mrk 501 and Mrk 
421, and 10^^ ergs s"^ in PKS 2155-304 |55]. These AGN are nearby BL Lac objects that 
correspond to FR I radio galaxies like Cen A seen along the jet axis, though Cen A may 
have the added advantage, in terms of UHECR production, of an external radiation field 
to enhance photohadronic processes in the jet [33, |38]. It will be interesting to apply 
this technique to a sample of radio galaxies, including Perseus A (3C 84, NGC 1275; 
z = 0.018 or dL = 76 Mpc) [561 [57] , Cyg A [3l] , etc. 
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4.2. UHECR Acceleration 

The particle energy density of a cold relativistic wind with apparent isotropic luminosity 
L and Lorentz factor T = — P"^ at radius R from the source is Up = L/{AtiR^(3V'^c). 

If a fraction is channeled into magnetic field B' in the fluid frame, then RB'T = 
2eBL/(3c, implying maximum particle energies -E^ax — QB'{R/ PT), so 



i^.H^ on ,/eBL/1046 ergs s-i 
Ema. ^ E'^^J ^ ) ^ 2xlO^°Z^^ '^j^ eV . (17) 

This simple, optimistic estimate begs the question how to transform directed particle 
kinetic energy into magnetic field energy in a cold wind. In other words, it is simply a 
dimensional analysis that contains no physical basis, whether for the source that emits 
a wind with such large apparent powers, or for the particle acceleration mechanism. 

Within the picture of particle acceleration through Fermi processes, a maximum 
particle energy in colliding shells can be derived in a straightforward manner. The 
underlying limitation of Fermi acceleration, whether first- or second-order, is that a 
particle cannot gain a significant fraction of its energy on a timescale shorter than the 
Larmor timescale. A colliding shell picture due to inhomogeneities in the relativistic 
wind realistically applies to the inner jets of radio galaxies and blazars, and to GRBs. 

Consider the ejection of two shells, with shell a ejected at stationary frame times 
< < At^,a, and shell b at times t*d < < t*d + /S.t^.b. The coasting Lorentz factor, 
wind luminosity (assumed constant during the duration of shell ejection), and energy 
are Ta{b),L^.^a{b) and £^*,a(fe), respectively, with stars referring to the stationary jet frame. 
For a colhsion, pr = '^a/'^b < 1- For relativistic winds, i.e., 1, the colhsion radius 

is rcou = 2F^c(t^,rf — At^a) when pr 1- When F^ > F^ ^ 1, the shocked fluid Lorentz 
factor F ^ 1. 

Collisions between relativistic shells divide into a number of cases depending on 
whether the forward shock (FS) and reverse shock (RS) are relativistic or nonrelativistic 
|58j . We consider the case of a nonrelativistic reverse shock (NRS) and relativistic 
forward shock (RRS), which is probably most favorable for particle acceleration. Thus 
the FS Lorentz factor F/ ^ 1. The magnetic field of the forward- shocked fluid in the 
shocked fluid (primed) frame is 



_£f^ PjMh, , (18) 

' r col J- a V C 

and e^j, the €3 parameter for the FS, is familiar from blast wave studies (e.g., [59|). 

The maximum particle energy in the comoving frame is -E^ax — ^^E'^cAt'^., where 
At'^ is the comoving duration when particles are undergoing acceleration. If this is 
equated with the time that it takes for the FS to pass through shell a, then 

At; = ^M^JfA , (19) 



From this, we obtain 



771 FT?' — 

^max — J- -^max — 
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Ze r 2tBfLa ^ „ . .„9n„ r 



-S 2-4 X Vf'Z ^,f^ eV. (20) 



(itt-0 

The factor F/F^ is at best of order unity, so this expression shows that for protons to 
reach ^ 10^° eV cosmic ray energies through Fermi processes, it is essential to consider 
sources with apparent isotropic luminosities ^ 10^^ ergs s~^. 

Acceleration of UHECRs in Cen A is therefore in principle possible during powerful 
episodes of jet activity. Acceleration of UHECRs can also take place in colliding winds 
of moderately powerful blazars with L ^ 10^^ ergs s~^ where strong Doppler coUimation 
takes place. 



4.3. Deflection of UHECRs 

The milliarcsec-scale radio jets in Cen A show moderate asymmetry, consistent with Cen 
A's radio jet being mildly relativistic and misaligned by ~ 60° j60j. Off-axis radiation 
beaming factors would conceal on-axis Cen A blazar-type flares with apparent powers 
L ^ 10^^ ergs s~^. These events would eject cosmic rays to ^ 10^° eV energies into the 
~ 100 kpc X 500 kpc radio lobe structure. If the equipartition field characterizes the 
large-scale magnetic field of the lobes (i.e., Ni^y ~ 1), then 60 EeV UHECR protons 
with tl = 65 kpc (see eq. P), could be deflected by the magnetic field of Cen A's lobe, 
or in the lobes of other radio galaxies or BL Lac objects, as suggested in arrival direction 
maps [28] . 

If UHECRs are formed through a neutron beam, then they will travel on average 
~ 500[-E'/(60 EeV)] kpc before decaying. Thus UHECRs will be deposited throughout 
the radio lobe, with some not decaying until outside the radio lobe structure. Thus 
both radio galaxies and blazars can be sources of UHECRs. For distant sources, 
though, energy losses will reduce the arrival energies of UHECRs, so that searches 
for enhancements of UHECRs towards specific sources should also take into account the 
EBL-dependent horizon energy, as illustrated in Fig. [2l 

The proton Larmor radius is formally tl = IIOE20 / BpQ Gpc (eq. [1]), and 
photomeson interactions with the CMBR and EBL take place on a horizon distance 
scale ~ rhor{E2o), where the received proton has lost ~ 63% of its energy. If the proton 
is hardly deflected out of the beam, so c/ ^ tl, then we should expect a pulse pile- up 
at a specific energy rhor{E2o) = dL. For PKS 2155-304 aX z = 0.116 or a propagation 
distance d = 400 Mpc (Fig. [2]), searches should be made for enhancements of UHECRs 
with EcR ^ 40 EeV, which is insensitive to the level of the EBL. If BL Lac objects 
accelerate UHECRs, then enhancements at Ecr ~ 40 - 60 EeV can be searched for 
from Mrk 421 and Mrk 501. Failure to detect a signal would constrain the minimum 
IGM field and number of inversions, or call into question a radio/7-ray galaxy origin of 
UHECRs. 

In Fig. [3] we plot the arrival directions of 27 cosmic rays from Auger (blue circles) 
and 13 cosmic rays from HiRes (red circles), all with > 56 EeV energy and with 1° 
angular uncertainty in their arrival directions (reflected in the radii of the circles; cf. 
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Figure 3. Hammer-AitofF projection in galactic coordinates of UHECR arrival 
directions and directions to nearby prominent AGN; the direction to the Galactic 
Center is at the left and right extremities of this plot. Thick blue, red, and orange and 
magenta circles correspond to UHECRs with > 56 EeV from Auger, > 56 EeV from 
HiRes, and > 40 EeV from AGASA, respectively. The radii of the circles reflect the 
angular errors in the reconstructed cosmic ray directions, which is 1° for both Auger 
and HiRes and 1.8° for AGASA. The centers of the concentric green, purple and pink 
circles correspond to the named AGN directions. The pink and purple circles represent 
angular deflections from the source AGN of the arriving cosmic rays with 40 EeV and 
20 EeV energies, respectively, in the /iG galactic disk magnetic field, using eq. 
The green circles, on the other hand, represents angular deflection in an assumed 0.1 
nG intergalactic magnetic field, assuming no magnetic-field reversals. The dashed and 
dot-dashed curves correspond to the equatorial and super-galactic planes. 

Ref. [61]). We also plot the arrival directions of 58 AGASA UHECRs, consisting of 
24 with E > 56 EeV (magenta circles) and 34 with 40 EeV < -E < 56 EeV (orange 
circles). Their arrival directions are displayed with an average angular resolution of 
1.8°, as previously reported by the AGASA experiment [62]. The number of AGASA 
UHECRs may suggest, given AGASA's exposure relative to HiRes and Auger, an energy 
calibration discrepancy between these experiments, which should be considered in more 
detailed studies but is ignored here. Note that the arrival directions of the cosmic rays 
from the prominent nearby AGNs (labelled) are deflected by the galactic (pink and 
purple circles) and inter-galactic (green circles) magnetic field. We used two energies, 
40 EeV (pink circles) and 20 EeV (purple circles), to calculate deflections in the fj,G 
galactic magnetic field according to eq. ([2]). The deflections from the AGN near the 
galactic plane (horizontal line) are larger because of the larger magnetic field. 

If Cyg A, BL Lac, 1959-1-65 or Mrk 501 are the possible sources of UHECR, then 
we should expect a large scattering, due to the galactic magnetic field, in the arrival 
directions of UHECRs below the GZK energy. For distant sources, e.g., 0152-1-17, the 
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deflection in the intergalactic magnetic field, plotted here for 0.1 nG field using green 
circles, may be larger. However, reversals of the field orientation over ~ Mpc scale 
may reduce the deflection in the intergalactic magnetic field from the values plotted 
in Fig. [3l These deflections in magnetic fields and related uncertainties can prevent 
positive identification of the UHECR sources below the GZK energy, although one 
should expect a clustering effect as evident from the lower energy AGASA data in the 
northern hemisphere. 



4.4- Nuclear 7 rays from Cen A 's Radio Lobes 



Typical ~ 1 fiG magnetic fields in Cen A's lobes carry an amount of energy = 
4 X lO^'^VyiiJ^Q ergs cm~'^, denoting the Cen A lobe volume Viobes = 
71 cm For magnetic fields near equipartition, only ^ few Myr are required for 



j^gss QYgs. Thus the Cen A jet 



a jet power of ~ 10^^ ergs s~^ and total energy of 
activity is likely to operate intermittently. 

Assuming at least equal energy in nonthermal protons and ions as nonthermal 
electrons, and using the equipartition assumption to normalize particle number and 
energy, we write the total cosmic ray spectrum over the Cen A lobes as 

WcR p-2 



Ncr{i) 



mp (ergs j 



1 _ 

-L Jmax 



7 



for 1 ^ 7 ~ '~1max 1 



(21) 



with 7maa: ~ 10^^ and total cosmic-ray energy Wcr = lO^^W^ss ergs. A 5-function 
approximation for the inelastic nuclear cross section, dapH-,-y/dE^ = 2aT,ox{'y)S{E^ — 
^Ep), where O"7r0x(7) — 27 In 7 + 58/^^7 — 41 mb is the tt" inclusive cross section [63] 
and 7-ray secondary fractional energy ^ ~ 5 - 10%, gives a lobe-integrated ^ 1 GeV 
uFi, spectrum 



2nd^ 



WcR{p-2)a^o{-f) , 7 = 



E^ 



(22) 



so 



4 X 10-i2^P-i(p-2) 
0.942-P 



(ergs cm' 

ne 

10-4 cm-3 



-2 -l^ 



58 



"3.5 



a^ox{E^/^mpC^ 
100 mb 



For ^ = 0.05 and p = 2.3, the nuclear 7-ray flux from the lobes of Cen A is 



2.4 X 10-^4^58^. 

"3.5 



-1-0.3 



10-4 cm-3 



a^ox(20^, 
100 mb 



ergs 
cm^ s 



(24) 



with E^ now in GeV. Based on lack of observed internal depolarization and measured 
soft X-ray flux, Ref. [17] argue that thermal particle target densities ~ 10"^ cm"^. 

The number of source counts with E^{GeV) > Ei detected with the Fermi 
Telescope for ^ = 0.05 and p = 2.3 is S{> Ei) ^ Xl^tfcev fE^dE^E^-'^AciE^) 
PS 0.7(X/0.2)At(yr)(n£/10~4 cm~^)W^^E~^'^ , where the effective area of the Fermi 

Telescope is Ag{E^) = 8500/y^ cm^ below 1 GeV and Ag{E^) = 8500 cm^ at 
higher energies, and X ^ 1/5 in the scanning mode. The number of background 
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counts from the angular extent, AQ = 12/57.3^, of the lobes of Cen A, using as 
background the diffuse extragalactic 7-ray background of Galactic background [61], is 
B{> E^) = 280{X/0.2)At{yT)E-^■^ . Unless > 10~^ cm'^, this emission signature 
would be too faint to be detected with the Fermi Gamma ray Space Telescope. Thus the 
cosmic-ray induced nuclear emission in the lobes of Cen A is probably entirely negligible. 
The much more optimistic estimates in Ref. [65] are due to the assumption of a very 
soft injection index p = 2.7 extending from GeV to ZeV energies, which requires a much 
larger cosmic-ray power. 



4-.5. Thomson- Scattered CMBR 

The same electrons that make the radio synchrotron radiation also upscatter photons 
of the surrounding radiation fields. Using a 5-function approximation for Thomson 
scattering of quasi-monochromatic photons with mean energy m^c^eo and energy density 
Uo gives the Thomson-scattered uFy flux 

fl = (^) f:::V^H{y; 0,1), y = ; (25) 

the value of the Heaviside function H{y; 0, 1) = 1 when < y < 1 and H{y] 0, 1) = 
otherwise. For the CMBR, mgC^eo = 1.24 x IQ-^meC^ = 0.63 meV = 153 GHz and 
f/o = 4 X 10-13(1 + z)^ ergs cm'^ (recall that Eb = B/Bcr). Thus 

f^=^-^f:::^Hie;0,eT). (26) 

For the CMBR, the mean dimensionless energy of a photon that is scattered by an 
electron radiating most of its synchrotron emission at epk is 

ex = ^^J^ , (27) 
implying emission at 

- 43^^ ( ) MeV . (28) 



B^G VlOO GHz. 

The WMAP observations at 22.5, 32.7, 40.4, 60.1, and 92.9 GHz used by Hardcastle et 
al. [U] allowed these authors to predict the Thomson-scattered flux. As can be seen 
from this approximation, this flux extends only to the low-energy end ( ^ 100 MeV) of 
the Large Area Telescope frequency range on the Fermi Gamma ray Space Telescope 

m. 



5. Cascade 7 rays from UHE Neutral Beams 

For flares in BL Lacs exceeding ^ 10^^ ergs s~^, and the much more energetic flares 
in FR II galaxies and flat-spectrum radio quasars reaching apparent 7-ray powers 

§ The Galactic diffuse 7-ray background has a softer spectrum than the extragalactic background, so 
is probably not significant at i?^ ^ 1 GeV at the w +20° galactic latitude of Cen A. 
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^ 10^^ ergs s~^, we describe the sequence of events starting to produce cascade radation 
following the acceleration of particles in black-hole jets [38j. A neutral beam is formed as 
a consequence of photomeson interactions of the accelearted UHECR protons with the 
synchrotron and scattered radiation field in the inner jet. Some ^ 10% of the hadronic 
energy can be transformed into an escaping neutron beam, with a few percent in a UHE 
7-ray beam and a few percent into neutrinos [371 [38] . Escaping neutrons with energy 
En travel i?„/10^° eV Mpc before decaying into protons and low-energy /5-decay leptons 
and neutrinos. 



Figure 4. Mean-free paths 
for attenuation of 7 rays by 
the EBL (A^^), for Compton- 
scattering energy loss of rela- 
tivistic electrons or positrons 
with the CMBR in the Thom- 
son (At) and Klein-Nishina 
(Akn) regimes, and neutron 
decay (A„eu)- The photopion 
energy-loss pathlength for pro- 
tons, from Fig. [1] are shown 
for comparison. The parti- 
cle or photon energy is de- 
noted £;(TeV). Short-dashed 
lines give length for a lepton 
to be deflected by 0.16'_i rad 
g g g g g in a magnetic field from 10~^ 

10" 10 10 10 ^q2o 10 -lO-i^G, as labeled. Neutron 

GeV E (TeV) PeV EeV eV ZeV decay length is given by A„e«. 

The 7 rays formed by j9, n -|- 7 vr" 27 processes in the inner jet and 
particle beam are attenuated by CMBR photons through 77 —>■ on length scales 
A^^^(kpc)= 2Ep,v/H0AEp,v) kpc for Ep,v > 1 and A^';(kpc)^ 4^E^ exp{l/ Ep^v) 
kpc for Epev ^ 1, valid until Epev < 0.1, when absorption on the extragalactic 
background light from stars and dust dominates. Leptons with 7 = 10^79 Compton- 
scatter the CMBR photons, losing energy on length scales At = 0. 75/79 kpc when 
79 < 1, and Akn = 2.l79/[ln(1.879) - 2] kpc when 79 > 10. 

If the distance d to the source is smaller than the particle Larmor radius and the 
correlation length (formally, Ninv = 1), the cascade pairs are deflected out of a beam 
with opening angle 6 = 0.16'_i radian when 6ri^ < At, implying that the collimation 
of the electromagnetic cascade is preserved until 79 = OA^B^u/O^i. Before the beam 
disperses, CMBR photons are scattered to energies e ~ 10"^7^, or E^ = lOOi?_ii/0_i 
TeV. Fig. m shows that the ratio ^ 0.1 is unique in that Compton-scattered 

CMBR from the cascade emerges from behind the EBL absorbing screen to form a hard 
7-ray component at 10 GeV - 10 TeV energies. If Bjgm ~ 10~^^ G, then the cascade 
emission from sources like lES 1101-232, if UHECR sources, could make an anomalous 
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Figure 5. Rectilinear 
cascade formed by 
photohadronic processes 
from UHECRs produced 
by Cen A. Lower-energy 
blue family of curves 
shows the synchrotron 
emission from succes- 
sive generations of the 
cascade, and higher- 
energy family shows 
various generations 
of the Compton/7-7 
cascade. One-year 
Fermi Gamma ray Space 
Telescope/GLAST and 
50 hr VERITAS sensi- 
tivities are shown for 
comparison. 

hard 7-ray component. 

The high-energy cascade flux induced by rectihnear motions of UHECRs from Cen 
A towards us is shown in Fig. [51 This calculation is a one-zone model, because it 
assumes that the properties of the IGM remain unchanged between Cen A and the 
Galaxy. Cosmic ray propagation through the extended magnetized environment in the 
cocoon and radio lobes of Cyg A [31] can also produce a 7-ray signature, but Cyg 
A is beyond the GZK horizon, so only lower energy UHECRs can reach us from this 
source. (For UHECRs accelerated in galaxy clusters |67] or other accelerators see, e.g., 
[68].) This calculation is normalized to 2 protons from Cen A with E > 60 EeV for the 
exposure and area of the Auger observatory (eq. [6]). 

The parameters of this calculation are d = 3.5 Mpc, B = 10^^ G, and particle 
injection index a = 2.2. The proton injection spectrum is exponentially cutoff at energy 
Emax = 2 X 10^'^ eV. The EBL used is from Ref. [69]. Twelve cascade cycles are sufficient 
to reach convergence, and this is the number used to calculate the total Cen A spectral 
energy distribution induced by UHECRs interacting with intervening radiation fields. 
The synchrotron cascade radiation is colored blue. For the assumed magnetic field, the 
synchrotron flux is produced by 7 ^ 10^^ electrons and remains coUimated. However, 
the Compton flux between ^ 1 - 100 TeV does not take into account the deflection of 
electrons. It should therefore be considered as an upper limit that could be reached if 
B < 10-12 G. 

The solid curve in Fig. [5] shows the multiwavelength uF^, flux, but is clearly far too 
faint to be seen with present technology. In this case, the proximity of Cen A works 
against it being a bright hadronic source: there is not enough pathlength to extract 
a significant fraction of the UHECR proton's energy. In future work (Atoyan et al. 
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2009, in preparation), details of cascades formed by UHE leptons and 7-rays formed 
by photohadronic processes will be presented. The prospects of detecting more distant, 
aligned blazar jets, if they are sources of UHECRs, could be more favorable if B ^ 10~^^ 
G. UHECR production in luminous flat spectrum radio quasars could make anomalous 
hard components at multi-GeV energies to be detected with the Fermi Gamma ray 
Space Telescope. Such emission signatures could also be seen in GRBs [70l 171] . 

6. Summary and Conclusions 

We have considered the implications of the discovery [H [2] by the Pierre Auger 
collaboration of clustering of UHECR arrival directions towards Cen A and AGN in 
the SGP. To guide out thinking, we have introduced a model-independent definition 
of the GZK horizon distance, shown in Fig. [2] (cf. EQ]). Assuming that black- hole 
jets energize and inflate the radio lobes of radio galaxies, a new method to calculate 
jet power based on the radio spectrum, size, and jet/counterjet ratio is applied to Cen 
A data. The average absolute jet power of Cen A is found from this technique to be 
^ 10^^ ergs s~^. The apparent isotropic power in a small beaming cone could easily, 
therefore, exceed 10^^ ergs s~^ during flaring intervals. 

An apparent jet power at least this great is needed to accelerate UHECRs through 
Fermi processes. The coUimated UHECRs, consisting mainly of neutron-decay protons, 
can be deflected by the ~ /iC fields in the lobes of radio galaxies (for acceleration in 
the lobes of radio galaxies see, e.g., [72]). This can give the appearance that UHECRs 
are emitted from Cen A's radio lobes [28], or allow UHECRs to originate from more 
distant, misaligned radio galaxies like Cyg A, which itself has a highly magnetized, 
~ 20 yuG cavity within several hundred kpc of the central engine [31]. The enhancement 
of UHECRs in the directions to specific radio jet sources at different redshifts would 
give valuable information about the IGM magnetic fields and intensity of the EBL. 

The inner jets of radio galaxies, including Cen A, can make an escaping neutral 
beam of UHECRs. Because Cen A is a nearby radio galaxy with its closer jet misaligned 
by ~ 35° - 70° to the line of sight, the UHECR flux received from it is small, and the 
cascade radiation flux is not detectable with current instrumentation. Neither is the 
secondary nuclear 7-ray emission from cosmic ray interactions with the lobe's thermal 
gas if its density is ~ 10^^ cm"'^ [47J, though the pair halo flux could be detectable 
[73] . The most prominent multi-MeV radiation signature is due to the CMB photons 
Compton-scattered by the radio-emitting electrons to soft 7-ray energies [66l HTj. Use 
of the transient event class in analysis of Cen A data [75j could help pull out the low- 
energy, 10 - 100 MeV signal, which would be important for normalizing its magnetic 
field and total jet power. 

For UHECR blazar sources pointed towards us, anomalous hard cascade 7 radiation 
spectra is potentially detectable with the Fermi Gamma ray Space Telescope or ground- 
based 7-ray telescopes, as we show in detail in future work. Detection of anomalous 
7-ray signatures in blazars (or GRBs) could reveal emissions from UHECR acceleration 



Ultra high energy cosmic rays from black hole jets of radio galaxies 



19 



in these sources. Lack of association of UHECR arrival directions with radio galaxies 
could require that various other classes of extragalactic bursting sources be admitted, 
including long duration GRBs [761 EZ], low luminosity GRBs [78], or magnetars. For 
instance, Ghisellini et al. [79] argue that the Auger UHECRs are not actually associated 
with Cen A, but with the background Centaurus cluster at ~ 40 - 50 Mpc. They argue 
that the rapid discharge of young, highly magnetized pulsars, which are the progenitors 
of magnetars, might have sufficient power to explain UHECR origin. Here the future 
pattern of UHECR arrival directions found with the Pierre Auger Observatory is of 
great interest, and whether a concentration builds toward the direction to Cen A. 

Spectral signatures associated with UHECR hadron acceleration in studies of radio 
galaxies and blazars with the Fermi Gamma ray Space Telescope and ground-based 
7-ray observatories can provide evidence for cosmic-ray particle acceleration in black- 
hole plasma jets. Together with IceCube or a northern hemisphere neutrino telescope, 
observations of PeV neutrino and MeV - GeV - TeV 7 rays can confirm whether black- 
hole jets in radio galaxies accelerate the UHECRs. 
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